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More than 635 000 positional observations of planets and spacecraft of different types, 
mostly radiotechnical ones (1961-2010), have been used for estimating possible changes of the 
gravitation constant, the solar mass, and semi-major axes of planets, as well as the value of 
the astronomical unit, related to them. The analysis of the observations has been performed 
on the basis of the EPM2010 ephemerides of IAA RAS in post-newtonian approximation. 
The EPM ephemerides are computed by numerical integration of the equations of motion of 
the nine major planets, the Sun, the Moon, asteroids and trans-neptunian objects. 

The estimation of change for the geliocentric gravitation constant GMq has been obtained 
{GM Q ) /GM Q = (-5.0 ± 4.1) • 10" 14 in year (3a). 
The positive century changes of semi-major axes dj/aj have been determined simultaneously 
for the planets Mercury, Venus, Mars, Jupiter, Saturn provided with high-accuracy sets of 
the observations, as expected if the geliocentric gravitation constant is decreasing in century 
wise. Perhaps, loss of the mass of the Sun M the produces change of GMq due to the solar 
radiation and the solar wind compensated partially by the matter dropping on the Sun. 

It was been found from the obtained GMq change and taking into account the maximal 
limits of the possible Mq change that the gravitation constant G/G falls within the interval 

-4.2 ■ 10~ 14 < G/G < +7.5 ■ 10" 14 in year 
with the 95% probability. 

The astronomical unit (au) is only connected with the geliocentric gravitation constant 
by its definition. The decrease of GM & obtained in this paper should correspond to the 
century decrease of au. However, it has been shown that the present accuracy level of 
observations does not permit to evaluate the au change. The attained possibility of fining 
the GMq change from high-accuracy observations points that fixing the connection between 
GMq and au at the certain time moment is desirable, as it is inconvenient highly to have 
the changing value of the astronomical unit. 
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INTRODUCTION 

The whole array of high-precision observations of the planets and the development of 
modern planetary ephemeris form prerequisites for the study of very subtle effects, in par- 
ticular, a change of the geliocentric gravitational constant GM Q with time. The question 
of variability and the possible rate of the change of gravitational constant G is regularly 
raised and considered in some cosmological theories (Uzan, 2003; 2009). The Sun's mass 
Mq, can not be absolutely constant too. On the one hand, it decreases due to continuous 
thermonuclear reactions and production of the radiant energy, with the matter carried away 
by the solar wind. On the other hand, there is a regular drop of interplanetary substances 
on the Sun, including dust, meteoroids, asteroids and comets. 

In the history of the Sun there have possibly been the periods of positive and negative 
changes of the solar mass. In the initial period, shortly after the formation of the protosun 
and the beginning of nuclear reactions in it, the mass of the central body probably increased 
due to the continuing compression to the center of the initial cloud. During the formation 
of the planetary system, until the interplanetary space was cleared, the mass of the matter 
falling on the Sun was higher than the mass of the Sun reduced due to the light and cor- 
puscular radiation. Now traces of the original, protosolar cloud are only likely to remain 
on the periphery of the Solar system beyond Neptune, the issue of the quantifying changes 
of the mass of the Sun still remains open because of the uncertainty of the overall balance 
including the mass loss due to radiation and the matter carried away by the solar wind, and 
the mass increase due to the matter falling on the Sun, in particular, comets whose collision 
with the photosphere of the Sun was repeatedly recorded by the SOHO space observatory 
(http:/ /ares.nrl.navy.mil/sungrazer/). Any estimation is difficult because of the complexity 
of reliable estimation of falling matter mass as well as the time- varied intensity and the angu- 
lar distribution of the solar wind in space. In this paper, we attempt to obtain experimental 
estimates of the change of the solar mass, namely, the geliocentric gravitational constant 
GM Q from the analysis of observational data of motion of planets and spacecraft. 
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The change of the value of the astronomical unit (au) is related to the change of the 
geliocentric gravitational constant. The astronomical unit is close in its magnitude to the 
average distance from the Earth to the Sun, but, by its definition it is only related to the 
heliocentric gravitational constant GM Q . From the obtained estimate GM Q , we can estimate 
the possible change in time of the au value. This estimate can be compared with the direct 
change value of au obtained from observations. 

EXPECTED EFFECTS OF THE CHANGE OF THE SUN'S MASS 

Estimations of the mass change of the Sun and its rate repeatedly are cited in the papers 
related to the exploration of solar physics, the solar wind and radiation (e.g. Sunyaev, 
1986; Livingston, 2000). The Sun's luminosity L & somewhat varies during the eleven year 
cycle and the ~ 27-day rotation around its axis; however, the fluctuations L Q do not exceed 
0.1 4 0.2% (Frohlich, Lean, 1998; 2004). If we take the average total solar luminosity to be 
L & = 3.846 • 10 33 erg/s and the mass of the Sun M & = 1.9891 • 10 33 g (Brun et al, 1998), 
then the decrease of the mass of the Sun due to radiation as a fraction of the solar mass is 
equal to M = -6.789 • 1O" 14 M per year. 

The mass carried away with the solar wind was also repeatedly evaluated. The basic 
composition of the solar wind is as follows: approximately 95% — protons, 4% — nuclei of 
the helium atoms, and less than 1% — nuclei of atoms of other elements (C, N, O, Ne, Mg, 
Ca, Si, Fe) (Brandt, 1973; Hundhauzen, 1976). The total number of particles flying away 
every second, is approximate by estimated as 1.3 • 10 36 (Kallenrode, 2004). The flow of the 
solar wind affects the activity of the Sun, coronal mass ejections. Typically, the average loss 
per year through the solar wind is estimated as 2-lO~ 14 M (Hundhauzen, 1976; Hundhausen, 
1997; Meyer- Vernet, 2007), that is, less than a third of the mass loss due to radiation. There 
are estimates (2 4 3) ■ 1O~ 14 M per year (Sunyaev, 1986; Carroll, Ostlie, 1996; Livingston, 
2000), where a value of 3 • 1O _14 M can be considered the upper limit of the mass carried 
away by the solar wind. The cumulative effect of the relative annual decrease of the mass of 
the Sun due to radiation and the solar wind can be restricted by the inequality 

- 9.8 • 10" 14 < M /M < -8.8 • 10~ 14 . (1) 

The reverse process occurs due to the fall of the dust, meteor, asteroid and comet sub- 
stance on the Sun. The dusty environment can not make a significant contribution to the 
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mass of the material fallen. According to the current data, the density of the interplanetary 
dust decreases with distance from the Sun, so that at the distance greater than 3 au there is 
not dust practically, with the 2/3ds of the interplanetary dust concentrating in particles of 
1CT 5 4- 10~ 3 g, and the size of dust particles being mostly 1 4- 10 /im (Mann et al, 2010). The 
total mass of the dust matter is estimated approximately as 10 19 4- 10 20 g (Sunyaev, 1986). 
Even with the assumption that all the mass reaches the Sun in several thousands years, the 
rate of the dust fall-out particles will be less than 1O~ 16 M per year. However, the dust 
particles smaller than 2 microns are swept by the solar pressure, while the ones greater than 
2 microns move towards the Sun. The most part of the approaching dust sublimates within 
0.1 au (~ 2Oi? ) and can not reach the surface of the Sun. It is also necessary to note, that 
a substantial part of the dust is carried away by the solar wind to the periphery of the Solar 
system (Mann et al., 2010). Therefore, the possible rate of the dust component falling out 
on the Sun is much smaller than (10~ 16 4- 10~ 17 ) ■ M per year. 

The larger particles, meteoroids and asteroids may fall on the Sun. Studies show that 
there is a constant migration of asteroids with an opportunity to complete the orbit evolution 
by the collision with the Sun (Farinella and et al, 1994; Gladman et al, 1997). The total 
number of small bodies is very large; the number of bodies larger than 1 km is about 
1 million. A significant part of asteroids move within a region close to the orbit planes of 
major planets of the Solar system, mostly situated in the belt between the orbits of Mars and 
Jupiter. The current estimates of the total mass of the asteroid belt give (13 ± 2) ■ 1O _1O M 
(Pitjeva, 2010b), i.e. less than 10~ 3 mass of the Earth. For the ring to be able to exist 
for tens or hundreds millions of years, the fraction of the outgoing annually material should 
be significantly smaller than 10~ 7 4- 10~ 8 of the total mass of the asteroid belt, in case the 
main asteroid belt is not replenished from the outside. With the outgoing material falling 
on the Sun not regularly, we find that the upper limit of possible mass of the Sun drop-down 
material from the main belt is less than (10~ 16 4- 10 -17 ) • M per year. Thus, we obtain a 
significantly lower value (two to three orders of magnitude) than the decrease of the solar 
mass by radiation and solar wind, so in the solar neighbourhood and in the field of the main 
asteroid belt there is no sufficient interplanetary matter migrating to the Sun to be compared 
with the decrease of the solar mass due to radiation. 

The mass of the matter that can come from distant regions of the solar system, mainly 
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in the form of comets (Bailey et al 1992), is more uncertain. There are trans- Neptunian 
areas — the Kuiper Belt, a cloud of the Hills, the Oort cloud. Currently, a large number 
of comets is detected in the immediate vicinity of the Sun (sungrazing comets) using the 



LASCO coronagraph (http:/ /lasco-w ww.nrl.navy.mil/[ ) installed at the SOHO solar space 
observatory (Marsden, 1989; 2005). Comets close to and often passing near the Sun are 
not long-living. They can disintegrate into fragments, or completely "fall apart". An ex- 
ample is a large family of Kreutz comets (Sekanina, Chodas, 2007). Some of them en- 
ter directly to the dense layers of the Sun. On the pictures of the SOHO observatory 



(http://sungrazer.nrl.navy.mil/index.php) the death of small fragments of comets in the so- 
lar photosphere is regularly recorded (kamikaze comet). The falling rocky bodies are more 
difficult to register, as during their approach to the Sun the glowing gas tail like that of 
icy objects and comets, does not form and the rocky bodies remain invisible. The overall 
contribution of the visible and invisible objects may be significant, although a reliable esti- 
mate of the total mass of substance reaching the Sun is extremely difficult. Nevertheless, the 
upper limits to the total mass can be specified. Statistics of comets discovered by the SOHO 
observatory gives about 500 comets for 30 -j- 35 months, on the average 170 -j- 200 comets 
per year. Many of them completely evaporate during their passage through the lower layers 
of the solar corona. There are registered repeatedly instances when the comets reached the 
photosphere. We assume that all the detected comets "vanished" and their mass increased 
the mass of the Sun, an overestimated value will result. The relatively small comets with 
the size from tens to hundreds meters are usually recorded, but for the upper estimate, we 
assume that the diameters of their nuclei (d com ) are several kilometers (as for the comets, 
which approach the spacecraft). If we take the average value of d com = 5 km, the density of 
3 g/cm 3 and double the result due to the missing and invisible falling objects, the annual 
upper bound is 

M , com /M < +3.2 ■ 10~ 14 . (2) 

This estimate is comparable to values obtained for the mass loss of the Sun due to radiation 
and the solar wind, although it seems to be overestimated, as the large falling objects have 
not been actually registered. This value can be considered the upper limit of the possible 
increase of the solar mass due to the material falling in the form of comets, meteors, asteroids 
and dust. 
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Now one can point to the two sides of the common interval, probably overestimated, which 
should contain the value of M /M . To obtain the lower limit, let us take the maximum 
loss estimate due to the solar wind, and at the same time, the zero drop of the material on 
the Sun. To find the upper limit, we use the maximum estimate ([2]) for the material falling 
on the Sun and the assumption that there is no mass loss due to the solar wind. Then, we 
obtain 

- 9.8 ■ HT 14 < M /M < -3.6 ■ 1(T 14 . (3) 

There are the restrictions to be kept in mind when trying to estimate the experimental 
changes of the Sun's mass. 

THE INFLUENCE OF THE CHANGE OF THE SOLAR MASS ON OR- 
BITAL ELEMENTS OF PLANETS 

The change of the M solar mass must lead to the appearance of variations of elements of 
the planet orbits, but a small and monotonic change affects only some certain elements. The 
effect is sought from extremely the small expected change (of the order or less than 1O _13 M 
per year), so it is sufficiently to consider the influence within the framework of the two-body 
problem (the Sun and the planet), as it is usually done, because the correction due to the 
influence of other bodies in the weak effect associated with the change of the central mass, is 
still several orders of magnitude smaller. The two-body problem with the variable mass has 
a long standing history and goes back to papers at the turn of the XIX and XX centuries 
by Gulden (1884), Meshchersky (1893), Stromgren (1903), Plummer (1906), etc. (see the 
detailed overview by Polyakhova, 1989). Variants of the isotropic mass variation in the two- 
body problem without the appearance of reactive forces, when no a particle is passed any 
pulse, were considered by MacMillan, Jeans, Armellini, Duboshin, Levi-Civita. A similar 
problem arises considering a possible change in time of the gravitational constant G(t) under 
the Dirac's hypothesis (Dirac, 1938) leads to the change of the mutual attraction forces and 
accelerations between bodies, and the analogous equations of the two-body problem. The 
analysis of equations for the central field of a variable mass body within the framework of 
General Relativity is given in the paper by Krasinsky and Brumberg (Krasinsky, Brumberg, 
2004). 

If we denote fi(t) the product of G(M & + m), where m - the planet mass, the vector 
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equation of the relative motion for a body with mass m is written 



(4) 



In general, we can assume that the gravitational constant G, incoming in fi(t), may depend 
on time. Since the central field remains while fx{t) (jlj) changes, then the area integral remains 
too, which is obtained immediately if the left and right side of (jl]) to vector multiply on r: 

d 

r x r = or — (r x r) = 0, 

then 

r x r = c. (5) 

The flat motion follows from the existence of the vector area integral ([5]). Energy for the 
case of the time-dependent value of /j,(t) is changing and is not more the integral of motion. 
Taking into account the monotony and smallness of the change of fi(t), it was shown (Jeans, 
1924), that the invariant holds 

nit) * a(t) = const, (6) 

where a - is the orbital semi-major planet axis. Sometimes this relation is called the adiabatic 
invariant of Poincare- Jeans, so as the first conclusion (J6j) was made by Poincare (Poincare, 
1911). The assumption, that fi(t) varies rather slowly, is essential for derivation of this result 
(Gelfgat, 1965). In our case, the expected value for Sun | fi Q (t) / fi Q (t) |~ 1CT 13 per year is 
much orders of magnitude smaller than it is required by the Gelfgat's constraints. 

Since the plane of motion is preserved, then the elements (i, Q), determining the position 
of the orbital plane, do not change. It should be considered the dependence of semi-axis 
a = a(t), the eccentricity e = e(t) and the argument of pericentre g = g{t) of the osculating 
elliptical orbits from the characteristics of the n{t) change for instantaneous values of coor- 
dinates, velocities, and a mass. From ([6]) we find a relationship between the changes of fi{i) 
and a{t): 

* = -1 (7) 
/i a 

The relationship obtained from the integral of areas for the osculating orbit is 

H(t) * a(t) * (1 - e 2 ) = c 2 , where c=|c|, (8) 

therefore, using fl6]), the e eccentricity of the osculate orbit remains constant under the given 
conditions e = const (Jeans, 1925). Investigation of the change of the g pericentre position 



for the small and monotonic //(i) change was made in the paper (Kevorkian, Cole, 1996), 
where it was shown that under accepted conditions of the smallness and monotony of fi, the 
g value does not have a secular trend, and there may be only small oscillations with the 
small amplitude of order (/i//i) 2 . 

The change in time of a = a(t), as leads to the change of the period of the body m 
that is leaving from the position for the case fi = const, and the increase of deviation will 
depend quadratically on the time interval. 

The possible fi & (t) = GM e change in the Solar system should be appeared in a system- 
atic, progressive, although very small deviation of the body position on the orbit (that is 
their longitude) and the change of the ai semi-major axes proportionally to the //©(t) change 
with opposite sign ([7]). The fact that the value //(£) = G(M Q + m) includes the m mass of a 
planet does not change the situation, since taking into account the mass m of a planet leads 
to the correction by several orders less. 

Thus, when the area integral remains and the attractive force from the main body de- 
crease/increase monotone, the second body is moving along the trajectory gradually reced- 
ing/approaching from/to the central body. The relative increase of the distance is equal to 
the relative decreasing of mass of the central body and vice versa. An orbit is transforming 
gradually remaining identical to itself, and is of a spiral form. The GM Q (t) change does 
not lead to secular trends of the eccentricity and the longitude of perihelion. The semi-axis 
is century changing a = a(t). Thus it is necessary to look for the effect caused by possi- 
ble change in time of the heliocentric gravitational constant in the corresponding secular 
variation of the semi-axes of the planet orbits. 

OBSERVATIONAL MATERIAL, REDUCTION OF OBSERVATIONS 

More than 635,000 positional observations of planets and spacecraft of various types 
(Table 1), mainly radiotechnical (1961-2010) have been used to construct high-precision 
ephemerides of planets and to determine the change of the heliocentric gravitational con- 
stant. The very accurate observations are required to find the very small effects, and it is 
most important and desirable to have observational data for planets, close to the Sun and 
having shorter periods, in the first place, the data for Mercury and Venus. Radiotechni- 
cal measurements which began in 1961 and are continuing with rising numbers since, first, 
yielded two new types of measurements in astrometry: the distance and the relative speed, 

8 



and secondly, the accuracy of the measurements became several orders of magnitude greater 
than the accuracy of the optical observations. 

Table 1. The observations used 



Observation type Time interval Observation number 



Optical 
Radiotechnical 



1913-2009 
1961-2010 



57768 
577763 



Total 



635531 



For this reason, the ephemerides of the inner planets provided by radiotechnical observations 
(mostly, data of time delays) are based fully on these data. At present, the radiolocation 
of planet surfaces is not carried out, but trajectory data of various spacecraft that orbiting 
around planets or passing near them are received regularly. Accuracy of observations of 
ranging has improved from 6 km to several meters for today's data of the spacecraft. It is 
necessary to say that the ephemerides of the outer planets so far, mainly, are based on optical 
measurements since 1913, when at the Naval Observatory of USA the improved micrometer 
was introduced, and observations become more accurate (0''5). 

Table 2. The distribution of optical observation and rms residuals in mas, 1913-2009 



Planet 


Observation number 


a 


Jupiter 


13038 


190 


Saturn 


16246 


150 


Uranus 


11672 


188 


Neptune 


11342 


177 


Pluto 


5470 


141 


Total 


57768 





However, until now a complete rotation period of Neptune and Pluto is not provided by 
the observations. In addition to optical observations of these planets, for the construction 
of ephemerides and estimation of their parameters the absolute observation satellites of the 
outer planets are used, as these observation are more precisely, and practically free from the 
phase effect hard taking into account, which is in observations of the planets themselves. 



Modern optical data are the CCD observations, their accuracy reaches 0''05. The number of 
used optical and radio observations, their planet distribution, as well as mean square error 
of residuals of observations are given in Table. 1-3. 

Table 3. Distribution, time interval, number and rms residuals for radiometric observa- 
tions 



l lane i 


uaia ^ype 


1 Hilt; 111 I t;I Vol 


kj usei vanon numuei 


<j 


Mercury 


t [m] 


1964-1997 


937 


575 




r [m] 


1973-2009 


5 


31.0 




a, S [mas] 


2008-2009 


6 


2.1 


Venus 


t [m] 


1961-1995 


2324 


584 




Magellan dr [mm/s] 


1992-1994 


195 


0.007 




MGN,VEX VLBI a, 5 [mas] 


1990-2007 


22 


3.0 




VEX r [m] 


2006-2009 


28163 


3.6 




Cassini r [m] 


1998-1999 


2 


2.4 




Cassini a, S [mas] 


1998-1999 


4 


105 


Mars 


r [m] 


1965-1995 


54851 


738 




Viking r [m] 


1976-1982 


1258 


8.8 




Viking dr [mm/s] 


1976-1978 


14978 


0.89 




Pathfinder r [m] 


1997 


90 


2.8 




Pathfinder dr [mm/s] 


1997 


7569 


0.09 




MGS t [ml 


1998-2006 


165562 


1.4 




Odyssey r [ml 


2002-2008 


293707 


1.2 




MRO r [m] 


2006-2008 


7775 


1.6 




spacecraft VLBI a, 5 [mas] 


1984-2010 


136 


0.6 


Jupiter 


r [m] 


1973-2000 


7 


13.8 




spacecraft a, 5 [mas] 


1973-2000 


16 


5.0 




spacecraft VLBI a, 5 [mas] 


1996-1997 


24 


9.5 


Saturn 


r [m] 


1979-2006 


34 


3.5 




a, S [mas] 


1979-2006 


92 


0.4 


Uranus 


Voyager-2 r [m] 


1986 


1 


7.4 




Voyager-2 a, 5 [mas] 


1986 


2 


11.0 


Neptune 


Voyager-2 r [m] 


1989 


1 


22.9 




Voyager-2 a, 5 [mas] 


1989 


2 


3.7 


Total 






577763 





The most accurate and long series of observations are available for Mars for which space- 
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craft and landers were launched repeatedly. Radiotechnical observations relating to Venus 
are much smaller, there are spacecraft Magellan and Venus Express. A situation for obser- 
vations of Mercury is much worse. The Messenger spacecraft (NASA) on the orbit around 
it have just appeared and we have not these data. There were only one-time conjunctions of 
Mariner-10 (1974-1975) and Messenger (2008-2009) spacecraft and ranging for the Mercury 
surface (1964-1997). Situation should be changed after the new data from the Messenger 
spacecraft and from the future spacecraft BepiColombo (ESA, launch 2014) will be avail- 
able. There are a number of radiotechnical observations for Jupiter and Saturn: the data 
for several spacecraft for Jupiter, and data of the Cassini spacecraft for Saturn. For Uranus 
and Neptune there are one 3-D points (a, 5, R), resulting from the conjunction of Voyager-2 
with these planets. The data were taken from the JPL database (http:/ssd.jpl.nasa.gov/iau- 
comm4/), created by Dr. Standish, and now supported by Dr. Folkner and the VEX 
data sent due to kindness of Dr. Fienga, as well as supplemented by rows of American 
and Russian radar observation of planets 1961-1995, taken from various sources. Russian 
radar observations of planets, along with their references are stored on the site of IAA RAS 
http: / /www.ipa.nw.ru/PAGE/DEPFUND /LEA /ENG /englea.htm[ The brief description of 
astrometric radio observations can be found in Table. 2 by Pitjeva (2005). 

EPM2010 PLANETARY EPHEMERIDES, DETERMINED PARAMETERS 

This work is based on the EPM2010 planetary ephemerides of IAA RAS. Numerical 
ephemerides of motion of the planets and the Moon (EPM — Ephemerides of Planets and 
the Moon) began to create in the seventieth years of last century under the leadership of 
G.A. Krasinsky for support of Russian space flights, and has successfully developed since 
then. The version of the EPM2004 ephemerides has used to release the Russian 'Astronom- 
ical Yearbook" described in (Pitjeva, 2005), the version of the EPM2008 ephemerides in the 
paper (Pitjeva, 2010a). The EPM2010 ephemerides were constructed using more than 635 
thousands of observations (1913-2010) of different types. EPM2010 differ from the previous 
versions by the improved dynamic model of motion of Solar system bodies, adding the pertur- 
bation from the ring of Trans- Neptunian objects (TNO), the new value of the Mercury mass, 
defined due to the three encounters of the Messenger spacecraft with Mercury, improvement 
of reductions of observations with the addition of the relativistic delay effect from Jupiter 
and Saturn, and the expanded database of observations, including radiotechnical (2008 - 
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2010) and CCD (2009) measurements. 

Ephemerides were constructed by the simultaneous numerical integration of equations of 
motion of all the major planets, the Sun, the Moon, the largest 301 asteroids, 21 TNO, the 
lunar libration, taking into account the perturbations from the oblateness of the Sun and the 
asteroid belt, lying in the ecliptic plane and consisting of the remaining smaller asteroids, 
as well as the ring of the TNO rest at the mean distance of 43 au. The equations of motion 
of the bodies were taken in the post-Newtonian approximation in the Schwarzschild field. 
Integration in the barycentric system of coordinates for the epoch J2000.0 performed by the 
Everhard method over 400 years (1800-2200) by the lunar and planetary integrator of the 
software package ERA-7 (Krasinsky, Vasilyev, 1997). The accuracy of numerical integration 
was verified by comparing the results of the forward and backward integrations over the 
century of the time interval. The errors were at least order of the magnitude smaller than 
the accuracy of observations. Thus, the accuracy of the ephemeris is determined mainly the 
accuracy of observations and their reductions. 

In the basic version of the improved EPM2010 planetary ephemeris about 260 parameters 
are determined: elements of the orbits of the planets and the 18 satellites of the outer planets; 
value of the astronomical unit; three angles of the orientation with respect to the ICRF frame; 
13 parameters of the rotation of Mars and the coordinates for the three martian landers; 
masses of 10 asteroids, the mean density for the three taxonomic classes of asteroids (C, S, 
M), the mass and radius of the asteroid belt, the mass of the TNO ring, the ratio of the mass 
of the Earth and the Moon; quadrupole moment of the Sun (J 2 ) and 23 parameters for the 
solar corona of different conjunctions of the planets with the Sun; eight coefficients of the 
topography of Mercury and the corrections to the level surface of Venus and Mars, constant 
shifts for the three series of planetary radar observations and for 7 spacecraft; 5 coefficients 
for the additional effect of the phase of the outer planets. 

The accuracy of EPM ephemerides was tested by comparison with the observations (all 
residuals do not superior to their a priori errors), as well as comparison with the DE421 
(JPL) independent ephemerides, those are in a good agreement (Pitjeva, 2010a). 

After constructing the EPM2010 ephemerides to all the observations the some other 
parameters can be to estimate: the changes of the GM Q heliocentric gravitational constant, 
the G gravitational constant, semi-axes of the planet orbits, and the astronomical unit. 
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OBTAINED ESTIMATES OF THE M ,G,ai CHANGES 

The main problem of this case consists in the smallness of the effects that need to be 
revealed. It was impossible to do this before an appearance in recent years a quite large 
number of high-precision observations, including data from spacecraft. Accuracy of deter- 
mination of the parameters increased significantly also due to extension of the time interval 
for which there are high-precision sets of planet observations. 

The parameters G and GM & were fitted by the least squares method simultaneously with 
all basic parameters of ephemerides, but each separately, i.e. they are considered in different 
solution versions. If GM & is found then it is taken into account that the accelerations 
between the Sun and other bodies change, and the mutual attractions between other pairs 
of bodies remain. This differs from the situation when we find the change of G and when 
the forces between all bodies vary accordingly. It should be noted that for the version of the 
G definition from the planet motions, the main contribution makes by the Sun, since the 
equations of the planet motions include products of the masses on the gravitational constant, 
among them the term for the Sun (GM ) is the main one of several orders of magnitude 
more than the others. Therefore, separating the change of G from the change M with the 
dominant term of the GM & is impossible. In this regard, it is more correctly (and reliably) 
to determine from planet motions the change of GM Q instead of G or M & separately. 

Table 4. The secular change values of the semi-major axes for the 6 planets provided 
with the high-accuracy observations 



Planet 


a/ a (century x ) 


Correlation coefficients 
between a and a 


Mercury 


(3.30 ± 5.95)-l(T 12 


56.5% 


Venus 


(3.74 ± 2.90)-l(r 12 


95.8% 


Earth 


(1.35 ± 0.32)-l(T 14 


0.6% 


Mars 


(2.35 ± 0.54)-l(T 14 


0.4% 


Jupiter 


(3.63± 2.24)-10" 9 


20.2% 


Saturn 


(9.44± 1.38)-KT 10 


35.9% 



Table. 4 shows the values obtained for the relative change of semi-major axes of the 
planet orbits. The most accurate results connect with availability of observations obtained 
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with using radiotechnical equipment, in particular using the spacecraft observations and the 
duration of larger time intervals. Accordingly, the most reliable relative values oi a/a have 
been received for them. These are the results for all the inner planets from Mercury to Mars. 
For Jupiter and Saturn the accuracy of a/ a is less. Values for the planets, not provided by 
radiotechnical data are unreliable. It is important that all the values obtained for the planets 
from Mercury to Saturn show the positive values of the a/ a ratio, i.e., indicate the decrease 
in time of the GM & heliocentric gravitational constant ([7]). 

The change of the GM & heliocentric gravitational constant has been determined from 
fitting all observations: 

(GM & )/GM Q = (-5.04 ± 4.14) • 1(T 14 per year (3a). (9) 

This was made similarly to find G/G: 

G/G = (-4.96 ± 4.14) • 1(T 14 per year (3a). (10) 

The closeness of the results ([9]) and ( ITOl) is not surprising, since while finding the G/G, 
when the forces between all pairs of bodies change, the effect of the central body is dominant 
and again the effect of GM G is found practically, instead of G. From the result ((9]) obtained 
for GM Q , it is possible to estimate the G value using the relation 

fi Q /fx Q = G/G + M Q /M Q . (11) 

This relation is valid with the 95% (2a) probability: 

- 7.8 • 10~ 14 < G/G + M & /M & < -2.3 • 10" 14 year" 1 . (12) 

Hence, using the limits (|3]) found for the value M /M , we obtain the G/G value with 
the 95% probability is within the interval 

- 4.2 • 10~ 14 < G/G < +7.5 • 10~ 14 per year. (13) 

Note, the G/G estimate, obtained in 2004 from to the lunar laser ranging (Williams et al, 
2004), which in any case is not complicated by the possible change of the solar mass, gives 
the following values of the limits for the gravitational constant change: G/G = (4±13)-10~ 13 
per year. The estimate of (GM Q )/GM Q obtained by us (jHJ) has the opposite sign and its 
value is an order of magnitude smaller. 
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The obtained change of GM Q , most probably, is related to the change of the Sun's mass 
Mq, rather than to the G change. Thus, we have 

M Q /M Q = (-5.0 ± 4.1) • 1(T 14 per year (3a). (14) 

Note that this value hits the limitation interval ([3]) for M & /M Q . The obtained change of 
GM Q OH]), possibly, reflects the balance between the mass loss due the radiation and the 
solar wind and the falling material contained in comets, rocky debris and asteroids, which 
do not produce the visible glowing gas tail. 

THE POSSIBLE CHANGE OF THE ASTRONOMICAL UNIT 

The change of the astronomical unit is connected with the change of the heliocentric 
gravitational constant. The astronomical unit, although is close in magnitude to the average 
distance of the Earth from the Sun, but by its definition (resolution MAS 1976 - IAU, 1976) 
is connected with the heliocentric gravitational constant: 

GM Q [m 3 s- 2 } = k 2 ■ au 3 [m 3 ] /86400 2 [s 2 ], (15) 

where k = 0.01720209895 is Gaussian gravitational constant. Currently, au is determined 
from ranging data with very high real accuracy, allowing to deduce the value of the he- 
liocentric gravitational constant from the of value au = (149597870700 ± 3) m, using 
the relation flT5j: GM Q = (1327124404 ± l)[£;m 3 s~ 2 ] (e.g., Pitjeva and Standish, 2009), 
which coincides with the value GM Q , proposed by W.Folkner and obtained by the same 
method. These values of au and GM & for the TDB time scale were approved at the 
XXVII IAU General Assembly (2009) as the best current values of astronomical constants 
( |http://maia.usno.navy.mil/NSFA2/NSFA_cbe.html[ ). 

In the paper by Krasinsky and Brumberg (Krasinsky, Brumberg, 2004) from raging 
data 1961 - 2003, using a numerical theory of planetary motion, about coinciding with the 
EPM2004 (Pitjeva, 2005), the authors obtained the secular increase of the astronomical 
unit au = 15 m per century, which should correspond to the increase of the heliocentric 
gravitational constant 

GM Q j GM Q ~ 3 • 10~ 12 per year. (16) 

The positive change of au should correspond to the decrease of semi-major axes of planet 
orbits, and not vice versa, as sometimes this is claimed, and alternative theories of gravity 
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(Miura T. et al, 2009; Nyambuya G.G., 2010) are even constructed on this incorrect basis. 
Such the large positive change ( JT6l) of the does not correspond to estimations of physical 
processes in the Solar system (the solar radiation and wind, the matter falling on the Sun), 
and also to the estimate (9) obtained in this study: GM Q /GM Q = —5.04 ■ 10~ 14 per year. 
However, authors considered themselves that the increase of au and the heliocentric gravita- 
tional constant (JT6l) are rather parameters of agreement than the real change of the physical 
parameters. 

Analysis of the obtained results based on the observations described in this paper, and 
the EPM2010 ephemerides, shows that the present level of observational accuracy does 
not permit to evaluate the au change. In the paper by Krasinsky and Brumberg the au 
change was determined simultaneously with all other parameters, specifically, with the orbital 
elements of planets and the value of the au astronomical unit itself. However, at present it is 
impossible to determine simultaneously two parameters: the value of the astronomical unit, 
and its change. In this case, the correlation between au and its change au reaches 98.1 %, 
and leads to incorrect values of both of these parameters, in particular, gives au the order 
of 15 m per century. 

Without the simultaneous determination of au and au, i.e. if only the change of the 
astronomical unit is estimated, together with other parameters, the au value is about 1 m 
per century, and does not exceed its formal uncertainty, thus it is determined: 

au = (1.2 ± 3.2) m/cy (3 a). 

Furthermore, including or excluding the au value from the number of the solution param- 
eters does not change the observation residuals, the mean error of the unit weight is also not 
changed (Act ~ 0.2%), so there is no reason to assume that au is the necessary parameter 
of agreement, and include it in the number parameters to be estimated. 

The modern accuracy has approached to the level when it is possible to estimate the 
change of the heliocentric gravitational constant GM & , therefore it is desirable to specify 
the definition of the astronomical unit, for example, by fixing the connection between GM & 
and au at the certain time moment, as it is inconvenient highly to have the changing value 
of the astronomical unit. 

CONCLUSION 

Modern radiotechnical observations of the planets and spacecraft having the meter ac- 
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curacy (relative error of 1CT 12 -j- 10~ n ) make it possible to obtain estimates of very small 
effects in the Solar system. The significant progress is related to several factors: increase in 
accuracy of reduction procedures for observations and in dynamical models of planet motion, 
as well as improvement of the quality of observational data, increasing their accuracy and 
the time interval in which these observations are obtained. 

The results obtained indicate on the decrease of the heliocentric gravitation constant per 
year at the level 

GM Q I GM Q = (-5.0 ± 4.1) • 1(T 14 (3a). 
For the gravitation constant it is found that value of G/G falls in the interval 

-4.2 • 1(T 14 < G/G < +7.5 • 1(T 14 

with the 95% probability. 

The obtained change of GM & seems to be due to the change of the solar mass M & , 
rather than the G change and reflects the balance between the loss of the solar mass due to 
by the radiation and the solar wind and matter falling on the Sun. It is possible to make the 
cautious conclusion that at present in the Solar system there is still the significant effect of 
matter falling on the Sun, that compensate partly the effect of reducing the solar mass due 
to the radiation and the solar wind. 

In the future, the connection between GM Q and au (15) should be fixed at the certain 
time moment, as it is inconvenient highly to have the changing value of the astronomical 
unit; moreover, it should be to define the changes GM Q , M , G rather than the change of 
the astronomical unit. 
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